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1 .  Background 


Since  1973  the  author  has  been  examining  the  statistical 
properties  of  Arctic  sea  ice  thermal  infrared  (TIR)  tempera¬ 
ture  measurements  obtained  with  VHRR  satellites  (LeSchack,  1974 
1975;  1976a),  The  thrust  of  that  research  was  to  determine 
whether  it  is  feasible  to  objectively  (i.e.,  numerically)  iden¬ 
tify  different  sea  ice  types  and  develop  a  system  for  automated 
mapping  of  sea  ice.  At  a  conference  where  the  results  of  the 
above  studies,  along  with  those  of  submarine  under-ice  analyses 
were  presented  (LeSchack,  1976b),  Buck  (1975)  made  the  observa¬ 
tion  that  it  would  be  valuable  to  his  under-ice  acoustic  work 
if  some  correlation  could  be  found  between  the  satellite  and 
submarine  data  sets,  both  recorded  in  the  eastern  Beaufort  Sea. 
Although  the  idea  of  making  a  large  area  correlation  between 
submarine  and  airborne  or  satellite  data  of  sea  ice  was  not 
new  (see  for  example,  Barnett,  1966),  this  was  the  first  time 
that  two  data  sets  usable  for  this  task  became  available.  Ac¬ 
cordingly,  with  both  the  availability  of  the  necessary  data 
and  the  impetus  provided  by  a  concrete  use  for  any  correla¬ 
tions  developed,  work  was  begun  in  1976  to  investigate  the  re¬ 
lationship  between  satellite  and  submarine  under-ice  roughness 
data . 


2 .  Under- Ice  Roughness 

The  variable  that  we  wish  to  ultimately  predict  is  under 
ice  roughness,  since  this  variable  strongly  influences  acous¬ 
tic  propagation  beneath  the  ice.  Several  mathematical  defini¬ 
tions  of  under  ice  roughness  can  be  formulated;  however,  only 
two  forms  appear  significant  to  our  problem.  These  are  the 
root-mean-square  (RMS)  under-ice  depth  and  the  standard  devia¬ 
tion,  a,  of  the  under-ice  depth,  respectively  defined  as: 


(1) 


(2) 


x  =  value  of  under  ice  depth  in  m  (depth  of  ice  from 
water  line) 

N  =  number  of  depth  values  along  a  profile 


At  the  outset  of  this  work,  RMS  under  ice  roughness  was  the 
parameter  of  choice.  This  was  because  not  only  is  it  a  measure 
of  roughness,  but  it  appears  to  be  a  good  indicator  of  overall 
ice  deformation  for  a  given  ice  surface  area,  since  any  signifi 
cant  departure  from  the  undeformed  equilibrium  ice  depth  over 
the  Arctic  Ocean  (about  3  m)  can  occur  only  through  building 
ice  ridges  and  keels  or  opening  of  leads  and  polvnyas. 

As  the  work  progressed,  however,  the  standard  deviation,  o, 
of  the  ice  depth  was  suggested  by  Buck  (1977)  as  being  more 
meaningful  to  the  acoustic  attenuation  problem  than  was  the  RMS 
ice  depth;  therefore,  in  subsequent  work,  o  was  used.  As  we 
shall  see,  there  is  a  close  correlation  between  these  measures 
of  under  ice  roughness,  as  well  as  the  "ice  roughness  factor" 
suggested  by  Buck  (1979). 

The  under-ice  data  used  for  the  early  portion  of  this  work 
were  collected  on  three  submarine  cruises  (Figure  1)  and  were 
discussed  in  detail  in  LeSchack  and  Chang  (1977).  The  SSN 
SARGO  collected  data  during  February  1960.  and  two  cruises 
were  made  by  the  SSN  SEADRAGON  during  the  summers  of  1960  and 
1962.  Selected  RMS  ice  depth  values  from  each  cruise  are 
plotted  on  Figure  1.  From  this  plot,  a  general  trend  of  in¬ 
creasing  RMS  ice  depth  from  West  to  East  appears,  with  maximum 
values  occurring  along  the  Canadian  Archipelago  (Figure  2). 
Comparison  of  these  data  distributions  with  the  distributions 
of  "ridging  intensity"  values  of  Hibler  et  al  (1974)  in  the 
Western  Arctic  Basin  shows  that  the  RMS  ice  depth  values  can 
be  well  delineated  by  the  three  ice  ridging  provinces  which 
they  call  Beauf ort-Chukchi ,  West  Central  Arctic  Basin,  and 
Archipelago.  Since  the  data  of  Hibler  et  al  (1974)  were  re¬ 
corded  a  decade  after  the  submarine  data  were  gathered,  a 
stable  pattern  of  ice  deformation  provinces  over  time  is  sug¬ 
gested.  Additional  under-ice  data  recorded  by  the  SSN  GURNARD 
in  1976  (to  be  discussed  later)  and  the  statistical  analysis 
of  the  1960-1962  data  (LeSchack  and  Chang  1977)  are  also  con¬ 
sistent  with  the  suggested  stable  under-ice  roughness  pattern 
of  the  Arctic  permanent  pack. 

As  mentioned  above,  preliminary  comparisons  between  stand¬ 
ard  deviation^  mean  ice  depth,  RMS  ice  depth  and  "ice  roughness 
factor"  have  been  made.  Figure  3  shows  the  interelationship  be¬ 
tween  the  mean  ice  depth,  the  RMS  ice  depth  and  the  standard 
deviation  for  the  data  recorded  by  the  SSN  GURNARD.  In  Figure 
4  the  standard  deviations  of  ice  depths  recorded  on  the  1960 
cruise  of  the  SSN  SARGO  have  been  compared  with  the  "ice  rough¬ 
ness  factor”  of  Buck  (1979),  This  factor  is  defined  as  the  pro 
duct  of  the  number  of  keels  equal  or  greater  than  2  m  per  50 
km  and  the  square  of  the  mean  ice  depth  for  the  same  50  km  ice 
profile  section.  This  factor  is  similar  to  the  "ice  ridging 
intensity"  of  Hibler  et  al  (1974). 
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FIGURE  1:  Cruise  tracks  of  the  SSN  SARGO  in  February  i960  and 
the  SSN  SEADRAGON  in  the  summers  of  1960  and  1962.  A  sampling 
of  RMS  ice  depths  (m)  along  each  track  is  included.  Each  data 
point  was  derived  from  a  segment  of  under-ice  profile  approxi¬ 
mately  8  nautical  miles  in  length.  (Modified  from  Lyon,  1963) 
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ICE  ROUGHNESS  FACTOR  (X10* 


FIGURE  4:  Plot  of  Ice  Roughness  Factor  (see  text)  vs  standard  deviation  about  mean  ice 
depth.  Each  data  point  is  derived  from  a  50-km  profile  section  from  the  1960 
winter  cruise  of  the  SSN  SARGO. 


As  can  be  seen  from  both  Figures  3  and  4,  there  appears  to 
be  a  close  relationship  between  all  measures  of  under  ice  rough¬ 
ness  that  have  been  used  by  various  researchers.  As  a  result, 
the  three  following  examples  of  correlations  of  TIR  data  with 
under-ice  roughness  have  been  presented  as  they  were  originally 
computed,  the  first  two  with  RMS  values,  and  the  last  with  values 
of  standard  deviation,  since  the  object  of  the  work  to  date  has 
only  been  to  establish  the  existence  of  a  functional  relationship 
between  TIR  and  under-ice  roughness  data. 

3.  Preliminary  Correlation  of  RMS  Ice  Depth  with  Satellite 

TIR  Temperature  Data _ 

Evidence  from  previous  research  (LeSchack,  1974,  1975,  1976a), 
revealed  that  when  a  statistically  sufficient  number  of  surface 
temperature  measurements  were  made  of  the  undeformed  multi-year 
ice  surface  in  the  Arctic,  using  airborne  or  satellite  TIR  scan¬ 
ners,  the  measurements  would  fall  into  a  normal  or  Gaussian  dis¬ 
tribution.  As  the  ice  deformed,  the  surface  temperature  dis¬ 
tribution  skewed  toward  the  side  of  the  curve  corresponding  to 
the  lower  temperatures.  The  temperature  distribution  appeared 
to  skew  more  with  the  greater  amounts  of  deformation. 

The  physical  basis  for  this  relationship,  to  be  discussed 
in  greater  detail  in  a  forthcoming  report,  appears  to  relate 
to  surface  ice  roughness  only  indirectly.  The  reasons  for  an 
increased  number  of  colder  infrared  temperature  measurements  with 
increased  ice  deformation  appear  to  relate  to  how  the  TIR  scan¬ 
ner  perceives  the  roughened  surface,  rather  than  to  any  direct 
sensing  of  lessened  heat  flow  due  to  thicker  ice  in  deformed 
areas.  From  an  analysis  of  the  sea  ice  surface,  the  following 
factors  seem  to  be  of  greatest  importance.  With  increasing  ice 
deformation : 

(1)  there  is  an  overall  increase  in  snow  catchment,  thus 
reducing  the  surface  temperature  directly  and  because 
snow's  emissivity  is  less  than  that  of  sea  ice,  re¬ 
ducing  also  the  TIR  scanner's  perception  of  the  sur¬ 
face  temperature. 

(2)  there  is  an  increase  of  shadowing  of  the  sun's  radi¬ 
ation  during  daylight  hours,  relative  to  an  undeformed 
surface.  These  shadow  areas  will  naturally  be  cooler 
than  irradiated  areas. 

(3)  there  is  an  increase  in  aerodynamic  surface  roughness. 
This  effect,  depending  on  conditions,  may  result  in 
either  slightly  increased  or  decreased  surface  temper¬ 
ature  . 

LeSchack  (1974,  1975)  have  shown  by  statistical  analysis 
of  radiant  temperature  data  distributions  for  given  areas  of 
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the  ice  surface,  as  obtained  by  the  NOAA  VHRR  satellite,  that 
an  objective  identification  can  be  made  of  the  various  surface 
types,  i.e.,  first-year  ice,  multi-year  ice,  water,  landfast 
ice,  land,  and  clouds.  LeSchack  (1975)  notes,  for  example,  that 
an  area  of  all  multi-year  ice  produces  a  normal  distribution  of 
temperature  data  points  with  a  mean  value  close  to  that  pre¬ 
dicted  by  the  model  of  Maykut  and  Untersteiner  (1971),  with  a 
skewness  of  zero  and  a  kurtosis  (peakedness)  of  about  3. 1  If 
increasing  areas  of  warmer  first-year  ice  are  included  in  the 
area  surveyed,  the  temeprature  distribution  skews  positively. 

If  the  area  includes  colder  clouds,  the  distribution  skews 
strongly  in  the  negative  direction  (<-0.20).  If  there  is  a  lot 
of  ridging  of  the  multi-year  ice,  the  temperature  appears  colder 
than  the  mean  multi-year  ice,  thus  producing  a  modest  negative 
skewness.  Therefore,  if  by  the  selection  techniques  discussed 
in  LeSchack  (1975)  a  la.rge  expanse  of  ice  surface  can  be  identi¬ 
fied  without  cloud  cover,  it  should  be  possible  to  predict  the 
general  extent  of  surface  ridging  by  analyzing  the  skewness  of 
the  temperature  distribution. 

Figure  5  is  an  IR  image  recorded  over  the  Western  Arctic 
Basin  on  17  April  1974.  The  area  enclosed  in  the  square  (ap¬ 
proximately  500  x  500  N.  miles)  includes  the  Beauf ort-Chukchi , 
West  Central  Arctic  Basin,  and  Archipelago  regions  outlined  by 
Hibler  et  al  (1974).  The  area  enclosed  on  the  image  has  been 
computer  processed  so  that  individual  arrays  of  radiant  temper¬ 
ature  data  can  be  examined  statistically.  Each  array  consists 
of  a  32  x  32  array  of  1-km  imagery  resolution  elements.  As  men¬ 
tioned  above,  it  has  been  possible  to  characterize  this  entire 
area  by  the  mean  temperature,  the  skewness  and  the  kurtosis  of 
each  32  x  32  km  array.  If  the  skewness  matrix  (Figure  6)  for 
the  area  shown  on  the  imagery  in  Figure  5  is  now  examined,  it 
can  be  seen  that  a  transect  can  be  picked  that  extends  more 
than  500  miles  in  length,  through  a  cloud-free  section  of  ice. 
This  transect  crosses  all  three  of  the  ridging  provinces  dis¬ 
cussed  by  Hibler  et  al  (1974)  and  furthermore,  comes  close  to 
paralleling  the  track  made  by  the  SSN  SARGO  cruising  beneath 
the  ice  in  February  1960. 

Figure  7  shows  a  section  of  the  SSN  SARGO  track  that  nearly 
coincides  with  the  transect  shown  in  Figure  6.  At  several  ran¬ 
dom  locations,  RMS  ice  depths  for  8-mile  segments  have  been 
computed  from  the  under-ice  profiles  recorded  by  the  submarine's 
upward-looking  sonar.  The  transect  shown  in  Figure  6  has  been 
divided  in  approximately  equal  area  segments  that  correspond 
closely  to  the  locations  for  five  of  the  RMS  ice  depth  values 
plotted  on  Figure  7  (the  value  associated  with  the  shear  zone 
segment,  i.e.,  the  large  shaded  area  on  Figure  6,  has  not  been 
used  since  this  area  contains  considerable  first-year  ice  with 


‘For  an  explanation  of  these  terms  see  Appendix  A. 
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FIGURE  5:  NOAA  VHRR  infrared  satellite  imagery  of  the  Western  Arctic  Basin 
recorded  on  17  April  1974.  This  area  has  been  numerically  analyzed  and  values 
of  the  skewness  of  the  surface  temperatures  are  presented  in  Figure  6.  All 
three  of  the  ridging  intensity  provinces  discussed  by  Hibler  et^  a^.  (1974) 
are  included  within  the  area,  which  is  approximately  500  x  500  n.  miles.  A 
NEE-SWW  transect  has  been  chosen  for  comparison  of  skewness  with  some  RMS 
ice  depth  values  (bold-face  numbers,  in  meters)  derived  from  the  February 


1960  under-ice  profile  data  recorded  by  the  SSN  SARGO  shown  in  Figure  7. 


WEST  CENTRAL  ARCTIC  BASIN’ 


MS--  -0.101 
RMS  =8.1 


MS=  -0.09 
RMS  =  7.2 


MS  =  -0.037 
RMS=  5.7 


MS  =  0.014 
RMS=  4.0 


MS=  -0.016 
RMS=  4.4 


FIGURE  6:  A  matrix  of  the  skewness  of  the  ice  temperature  distribution 
for  the  imagery  shown  in  Figure  5.  Each  skewness  value  was  derived  from 
a  32  32  km  array  of  temperature  data.  The  width  and  sections  of  the 

transect  were  arbitrarily  chosen  such  that  an  adequate  sample  of  skew¬ 
ness  data  could  be  compared  with  an  RMS  ice  depth  value  derived  from  a 
region  of  close  proximity.  The  mean  skewness  (MS)  and  RMS  ice  depth 
(RMS,  in  meters)  are  tabulated  for  each  section  at  right.  Banks  Island 
is  at  upper  righthand  corner.  Ridging  provinces  are  delineated.  The 
shaded  areas  have  a  skewness  >  0.50,  which  implies  much  first-year  ice. 
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resulting  high  positive  skewness  values).  When  the  data  pairs 
(i.e.,  the  skewness  and  the  RMS  values)  for  each  area  are  plot¬ 
ted,  a  clear  negative  linear  relationship  emerges  (Figure  8). 
For  this  to  be  a  meaningful  correlation,  we  must  assume  that 
the  ice  provinces  in  1960  were  essentially  the  same  as  they 
were  in  1974.  However,  the  work  by  Hibler  et  al  (1974),  the 
statistical  analysis  of  LeSchack  and  Chang  (1977),  and  the  con¬ 
sistency  of  the  SSN  GURNARD  data  of  1976  with  the  1960-62  data 
shown  in  Figure  2  all  suggest  that  the  type  of  ice  roughness 
stability  as  expressed  by  the  statistical  parameters  used  above 
is  possible. 

4.  Correlating  SSN  GURNARD  Data  with  Satellite  TIR  Data 


4.1  Introduction 


The  correlation  shown  in  Section  3  above  is  startling. 
However,  since  only  five  data  sets  were  available  for  the  re¬ 
gression  line  shown  in  Figure  8,  and  because  of  all  the  as¬ 
sumptions  as  to  ice  province  stability  that  are  necessary  be¬ 
fore  this  Figure  is  meaningful,  an  attempt  was  made  to  try  this 
correlation  again  with  a  new  set  of  data. 

In  April  1976,  the  SSN  GURNARD  gathered  under-ice 
roughness  data  along  three  straight-line  tracks,  each  about  200 
N.  miles  long,  in  the  vicinity  of  the  AIDJEX  Project.  The  cruise 
tracks  are  shown  in  Figure  2.  A  cloud-free  VHRR  IR  image  was 
obtained  for  this  area  of  the  Beaufort  Sea  for  8  March  1976,  about 
a  month  earlier  than  the  cruise.  The  cruise  tracks  were  laid 
out  on  the  rectified  imagery  (Figure  S)  and  16  x  16  km  squares 
of  TIR  data  (each  containing  256  pixcels)  along  these  cruise 
tracks  were  then  selected  for  analysis.  A  total  of  58  16  x  16  km 
(10  x  10  N.  mile)  data  squares  were  found  to  lay  along  the  cruise 
tracks.  A  statistical  analysis  of  the  temperature  distribution 
of  each  of  these  squares  was  then  made.  By  photo-interpretive 
techniques  a  determination  was  made  as  to  which  of  the  58  blocks 
were  all  multi-year  ice.  Some  22  blocks  were  chosen.1  At 
this  stage  of  the  research  it  appears  necessary  to  use  only  multi¬ 
year  ice  for  the  correlation  analysis;  mixtures  with  first-year 
ice  have  distorted  skewness  patterns  that  are  not  as  easy  to  cor¬ 
relate  with  under-ice  data  as  are  the  pure  multi-year  ice  dis¬ 
tributions  . 

Although  at  first  glance,  the  requirement  to  choose  only 
multi-year  ice  distributions  for  analysis  purposes  may  seem  limit¬ 
ing  in  the  context  of  conducting  a  large-scale  survey, in  fact, 
the  selection  process  appears  easy  to  do  by  computer,  an  import¬ 
ant  criterion  for  any  operational  use  of  this  prediction  tech¬ 
nique.  Figure  10  is  a  plot  of  TIR  temperature  skewness  vs  kur- 
tosis  for  the  58  16  x  16  km  data  blocks  chosen.  It  is  signifi- 


lAppendix  B  contains  a  tabulation  of  these  data. 
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SKEWNESS  OF  ICE  TEMP.  DISTRIBUTION  (  APRIL  1974  ) 


INI 


- - - - - - — - - . - 1 - 1 - 

4-°  5.0  6.0  7.0  8.0 

RMS  ICE  DEPTH  (M)  FEB  )  960 


FIGURE  8:  A  strong  negative  correlation  between  RMS  ice  depth  and  skew¬ 
ness  of  the  ice  surface  temperature  as  measured  by  NOAA  VHRR  satellites 
is  observed  for  the  few  locations  where  these  data  are  available. 
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cant  that  essentially  all  of  the  data  points  identified  by  photo¬ 
interpretation  as  being  pure  multi-year  ice  are  clustered  toget¬ 
her,  whereas  data  representing  a  variety  of  mixtures  of  multi¬ 
year  ice  and  first-year  ice  are  scattered  about  the  plot  in  no 
order  as  yet  discernible.  This  appears  significant  since,  as 
suggested  by  Figure  10,  an  algorithm  can  be  easily  constructed 
to  select,  with  high  probability,  all  the  pure  multi-year  ice 
data  in  a  given  population. 

4 . 2  Preparing  the  Two  Data  Sets 

The  most  obvious  criticism  of  the  correlation  presented 
in  Section  3  is  that  the  satellite  data  were  recorded  14  years 
after  the  submarine  data.  In  the  case  of  the  SSN  GURNARD  data, 
however,  satellite  data  of  the  area  were  available  for  the  pre¬ 
vious  month.  From  past  experience,  it  is  reasonable  to  assume 
that  during  the  March-April  period  there  will  be  little  physi¬ 
cal  change  in  the  ice,  although  as  can  be  seen  from  Figure  11 
a  westerly  driTt  might  be  expected.  Accordingly,  a  comparison 
of  satellite  and  submarine  data  taken  only  a  month  apart  can, 
with  adjustment  for  ice  movement ,  be  considered  to  be  made  be¬ 
tween  simultaneously  recorded  data. 

4 . 3  Registration  of  Satellite  and  Submarine  Data 

In  the  data  compared  in  Section  3,  a  broad  ice  surface 
(some  160  km  on  a  side)  was  analyzed  to  obtain  the  five  individ¬ 
ual  skewness  values  used  for  the  correlation  shown  in  Figure  8. 
Each  of  the  five  individual  skewness  values  represented  ice  con¬ 
ditions  for  a  whole  province  rather  than  a  small  area,  and  as 
a  result,  exact  registration  of  the  satellite  data  and  the  under¬ 
ice  data  with  which  it  was  compared  appeared  unnecessary.  In  the 
case  of  the  SSN  GURNARD  data,  however,  a  1:1  correlation  with 
the  satellite  data,  broken  up  into  16  x  16  km  elements,  is  being 
attempted.  Accordingly,  the  three  legs  of  the  cruise  have  been 
divided  into  adjacent  16  km  (10  N.  mile)  profile  segments  (Fig¬ 
ure  12).  The  mean  RMS  ice  depth  for  each  segment  was  computed 
and  tabulated  with  the  22  temperature  skewness  values  determined 
to  be  representative  of  pure  multi-year  ice  (Table  1).  The  skew¬ 
ness  data  have  been  plotted  in  the  geographical  position  along 
the  submarine  tracks  corresponding  to  a  direct  overlay  of  data 
sets;  i.e.,  the  16  x  16  km  satellite  data  blocks  have  nominally 
the  same  coordinate  position  as  the  submarine  data  beneath  them, 
notwithstanding  the  fact  that  some  ice  movement  is  anticipated. 

Examination  of  the  ice  drift  records  for  this  period 
shows  a  movement  of  10-20  N.  miles  in  the  western  quadrant 
(Thorndike  and  Cheung,  1977).  Accordingly,  it  was  felt  that 
some  shifting  of  the  data  sets  might  be  necessary  to  effect  data 
registration.  Figure  12  illustrates  the  data  matrix  configura¬ 
tion.  If  the  satellite  data  matrix  is  moved  north  one  block 
(16  km)  (and  we  assume  no  horizontal  deformation  of  the  ice  sur- 
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FIGURE  10:  A  plot  of  skewness  VS  kurtosis  of  temperature  distributions  derived  from  the  58 
16  x  16  km  data  squares  along  the  track  of  the  SSN  GURNARD  shown  in  Figure  9. 
(See  Appendix  B  for  data) .  Multi-year/first  year  ice  mixtures  have  been  de¬ 
termined  by  photo-interpretation  of  the  VHRR  image.  Note  that  all  of  the  pure 
multi-year  ice  data  cluster  at  the  bottom  of  the  plot 


FIGURE  11:  The  major  drifts  of  polar  ice  (after  Dunbar  and 
Wittman,  1963) .  Note  that  the  SSN  GURNARD  track 
falls  within  the  westerly-moving  ice. 
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SSN  GURNARD  Track 

16  x  16  km  squares  of  VHRR  TIR  data 

A  nominal  overlay  of  satellite  VHRR  data  on  the  RMS  under-ice  data 
Each  track  was  appropriately  divided  into  sections  of  sufficient 
length  to  match  the  corresponding  VHRR  data  blocks.  The  two  data 
matrices  can  be  shifted  until  the  best  registration  is  obtained. 


Table  1:  Comparison  of  Under  Ice  Roughness  with  Nominally  Overlain 
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face),  the  result  will  be  a  1  block  north  shift  along  track  1, 
approximately  a  1  block  northwest  shift  along  track  2.  and  no 
shift  along  track  3.  Similarly,  shifting  of  data  along  the 
three  tracks  resulting  from  moving  the  satellite  data  matrix 
to  the  NE ,  E,  SE,  S,  SW ,  W,  and  NW  can  likewise  be  accomplish¬ 
ed,  and  by  referring  to  Table  1,  a  different  set  of  data  pairs 
for  each  shift  can  be  derived  for  computing  correlation  func¬ 
tions.  Since  Thorndike  and  Cheung  (1977)  indicated  movement  of 
only  10-20  N.  miles,  a  maximum  matrix  shift  of  only  one  data 
block  (10  N.  miles)  was  considered  necessary  for  determining 
data  set  registration. 

4 . 4  Correlation  Functions 

Since  the  major  thrust  of  the  work  with  the  SSN  GURNARD 
data  was  to  verify  whether  the  correlation  presented  in  Figure  8 
is  real,  the  next  step  is  to  attempt  to  correlate  the  April  1976 
under-ice  data  with  the  March  1976  satellite  data  of  the  same 
area.  Using  the  two  data  matrices  presented  in  Table  1  and  al¬ 
lowing  either  no  ice  shift  or  a  10-N.  mile  shift  along  the  eight 
compass  points,  9  different  sets  of  22  data  pairs  can  be  gener¬ 
ated.  Although  in  the  case  of  Figure  8  the  small  number  of  data 
points  and  their  obvious  linear  relationship  makes  it  reasonable 
to  attempt  only  a  linear  fit  to  the  data,  a  linear,  parabolic, 
exponential,  logarithmic  and  power  fit  were  attempted  with  the 
SSN  GURNARD  data.  The  parabolic  fit  in  all  cases  provided  the 
best  correlation,  with  the  linear  fit  next  best.  The  9  sets  of 
correlation  results  are  presented  in  Table  2,  first  using  all 
22  data  pairs  and  then  eliminating  two  data  pairs  because  their 
skewness  values  appeared  too  high  to  be  truly  representative  of 
pure  multi-year  ice  (the  removal  of  these  two  pairs  is  based 
solely  on  an  intuitive  feel  for  the  data,  but  as  will  be  seen 
shortly,  appears  reasonable). 

Examination  of  the  correlations  in  Table  2  indicates 
that  only  when  the  satellite  data  are  shifted  SW,  W  or  NW, 
relative  to  the  submarine  data  are  there  any  significant  cor¬ 
relations.  In  Figure  13,  these  correlations  are  shown  graph¬ 
ically  when  the  parabolic  correlation  coefficients,  r, 
for  the  20  data  pairs  are  plotted  around  the  points  of  the 
compass.  Only  does  movement  in  the  western  quadrant  appear 
statistically  meaningful  while  movement  in  the  quadrants  away 
from  the  west  produces  monotonically  decreasing  correlation. 

This  is  consistent  with  the  westerly  ice  movement  in  the  March- 
April  1976  period,  i.e.,  the  satellite  data  matrix  had  to  be 
shifted  to  the  west  to  come  into  register  with  the  submarine 
data.  This  strongly  suggests  that  not  only  is  there  a  good 
correlation  between  the  satellite  TIR  data  and  the  submarine 
under-ice  roughness  data  but  that  the  correlation  is  meaning¬ 
ful  down  to  resolution  elements  of  16  km  (10  N.  mile)  on  a 
side . 
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TABLE  2;  Linear  and  Parabolic  Correlations  of  the 
Satellite  Data  Shifted  10  N.  Miles  with  Respect  to 
the  Submarine  Data 


22  Data  Points  20  Data  Points 

SHIFT _  Linear  (r)  Parabolic  jr) _ Linear  (r)  Parabolic  (r) 


None 

-.28 

.  30 

-.34 

.  36 

S 

-.40 

.41 

-.38 

.  39 

SW 

-.55 

.75 

-.67 

.89 

W 

-.65 

.81 

- .  66 

.88 

NW 

-.64 

.84 

-.65 

,89 

N 

-.41 

.42 

-.38 

.40 

NE 

-.18 

.20 

-.23 

,30 

E 

-.23 

.26 

-.19 

.  22 

SE 

-.28 

.33 

-.23 

,28 

.40 

N 


.39 


FIGURE  13:  A  strong  correlation  is  observed  in  the  western  quadrant  when 
the  correlation  coefficient  (r)  is  plotted  around  the  points 
of  the  compass  (See  Table  2) . 
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It  is  instructive  to  examine  a  plot  of  the  skewness 
RMS  data  for  one  of  the  shifts  that  produces  a  good  cor¬ 
relation.  Figure  1-1  is  a  plot  of  the  SW  shifted  data.  The 
two  points  that  have  been  eliminated  from  one  of  the  sets  of 
correlations  have  been  indicated  on  the  plot.  They  have  been 
eliminated  both  on  the  grounds  that  they  do  not  seem  to  clus¬ 
ter  vith  the  remaining  data  and  that  a  skewness  of  greater  than 
+0.3  probably  indicates  warm  clouds  over  the  ice  or  first-year 
ice  mixed  with  the  multi-year  ice.  either  of  which  adversely 
affects  the  relationship  of  skewness  to  ice  deformation.  The 
remaining  data  points  appear  generally  to  be  scattered  equally 
around  the  curve  fit  to  the  data. 

5.  Correlation  oi  Simultaneous  Airborne  TIR  Data  with 

Under-Ice  Data  by  HMS  SOVEREIGN _ 

5 . 1  Introduction 

A  third  study  was  then  conducted  that  shows,  perhaps 
more  graphically  than  the  others,  the  correlation  of  under-ice 
data  recorded  by  the  British  nuclear  submarine  HMS  SOVEREIGN 
between  18-21  October  1976  with  airborne  TIR  data  recorded  dur¬ 
ing  the  same  period  over  the  submarine  track  by  a  Canadian 
Forces  Argus  aircraft.  In  many  ways,  this  example  is  better 
than  the  previous  two.  The  airborne  TIR  data  are  recorded  in 
the  form  of  imagery  on  a  scale  large  enough  to  see  the  in¬ 
dividual  features  of  the  ice  surface.  Every  attempt  was  made 
to  cover  the  exact  track  of  the  submarine  within  hours  of  its 
passing  beneath  the  ice. 

Figure  15  shows  the  submarine  track.  The  data  dis¬ 
cussed  in  this  example  were  collected  along  Segment  AB .  Wad- 
hams  (1978)  divided  Segment  AB  into  10  sections,  each  approxi¬ 
mately  97  km  in  length.  They  are  numbered  as  in  Figure  15. 

The  airborne  TIR  data  used  in  this  correlation  were  recorded 
from  an  altitude  of  4000  ft  (1200  m)  between  1910  hrs  and 
2109  hrs  Z  on  19  October  1976  along  Segment  AB .  The  exercise 
was  referred  to  as  "Brisk  Laser  #1"  by  the  Canadian  Defense 
Research  Establishment. 

5 . 2  Data  Analysis 

5.2.1  Under-Ice  Data 

As  in  the  previous  Beaufort  Sea  examples  cited, 
the  object  of  this  work  is  to  derive  an  under-ice  roughness 
parameter  from  the  profile  data,  and  to  correlate  it  with  a 
parameter  easily  derived  from  the  airborne  TIR  data  that  physi¬ 
cally  relates  to  the  topside  roughness  of  the  ice.  Although  in 
the  previous  examples  the  RMS  ice  depth  was  used  as  the  under¬ 
ice  roughness  parameter,  recent  work  by  Buck  (1979)  indicates 
standard  deviation  (a)  from  the  mean  under-ice  depth  is  more 
meaningful  for  Arctic  acoustics.  Therefore,  a  will  be  the 
under-ice  roughness  parameter  used  in  this  example. 
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'IGURE  15:  Cruise  track  of  HMS  SOVEREIGN,  October  1976 


24 


Wadhams  (1978)  prepared  histograms  (probability 
density  functions)  of  the  under-ice  depths  recorded  in  each 
of  the  10  97-km  zones  along  Segment  AB.  These  are  shown  in 
Figure  16.  For  each  of  the  10  zones  P(H),  the  frequency  of 
occurrence,  is  plotted  against  the  ice  depth.  In  this  ex¬ 
ample,  the  standard  deviation  of  the  mean  ice  depth  was  com¬ 
puted  from  the  Wadhams ‘  histograms  using  the  following  re- 
lat ionship : 

fi  <di  -  a)2  +  f2  (d2  -  dr  +  .  +  fn  (dn  -  d)-- 


(f. 


-  l) 


where 


The 


c  =  the  standard  deviation 

f^  =  the  frequency  of  occurrence  for  each 

d^  =  the  ice  depth  at  the  center  point  of 


( i . e. ,  dx 
computed  values 


=  0.5  m,  d^  =  1.5  m,  d^  = 
for  a  are  listed  in  Table 


1  m  depth  cell 
each  1  m  cell 
2.5m,  etc . ) 

3  . 


5.2.2  Airborne  TIR  Data 

As  shown  by  LeSchack  (1974,  1975,  1976a),  when 
many  TIR  temperature  measurements  ai’e  made  over  a  sea  ice  sur¬ 
face,  each  ice  type,  i.e.,  multi-year  ice,  first-year  ice  and 
thin  ice,  has  associated  with  it  a  normal  or  Gaussian  distribu¬ 
tion  of  temperatures.  For  example,  if  a  given  ice  surface  is 
comprised  wholly  of  multi-year  ice.  an  adequate  series  of  sur¬ 
face  temperature  measurements  produce  the  familiar  "bell-shaped" 
distribution  curve.  If  there  is  multi-year  and  first-year  ice 
in  the  surface  area  measured,  two  slightly  overlapping  normal 
distribution  curves  are  produced.  When  thin  ice  is  also  in¬ 
cluded,  three  overlapping  curves  can  be  seen.  The  area  under 
each  curve  is  proportional  to  the  amount  of  surface  represent¬ 
ed  by  each  type. 


In  this  example,  the  airborne  TIR  data  were  re¬ 
corded  over  the  submarine  track  in  the  form  of  imagery  on  70  mm 
film,  The  imagery  used  was  recorded  between  1912  hrs  and  2040 
hrs  on  19  October  1976,  This  section  of  imagery  had  been  care¬ 
fully  annotated  by  Dr.  Bev  Young  of  the  Canadian  Department 
of  National  Defense  and  the  geographical  position  of  the  im¬ 
agery  can  be  easily  ascertained.  Between  1912  and  2040  hrs, 
aircraft  altitude  was  constant,  there  was  no  cloud  cover,  and 
the  IR  instrument  operator  left  the  gain  control  constant 
so  that  the  photographic  gray  scale,  although  not  temperature- 
calibrated,  maintained  the  same  temperature/gray  scale  relation- 
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FIGURE  16:  Probability  density  functions  of  ice  draft. 
Plots  are  numbered  according  to  divisions  seen  along 
Segment  AB  in  Figure  15.  (after  Wadhams , 1978) 
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ship.  During  this  period,  imagery  covering  sections  -1-8  of  the 
submarine  track  was  obtained. 

Five  sections  ol  imagery,  representing  the  ice 
surfaces  close  to  the  center  of  each  segment,  were  scanned  with 
a  Photometric  Data  Systems  Inc.  Series  1000  Microdensitometer1 
to  obtain  quantitative  values  proportional  to  the  ice  surface 
temperature.  Each  area  scanned  was  chosen  so  that  a  majority, 
if  not  all,  of  the  surface  was  multi-year  ice.  The  decision 
was  made  by  photo- interpret ive  means.  A  straight  scan  line  19 
cm  long  through  the  center  of  the  image,  parallel  to  the  line 
of  flight,  was  made.  A  circular  aperture  86  urn  in  diameter 
was  used  to  measure  density.  A  density  value  was  recorded 
every  86  um  along  the  scan  for  a  total  of  2210  samples  per 
scan . 


In  the  imagery  used,  the  lower  the  film  den¬ 
sity,  the  warmer  the  surface;  hence,  the  thinner  the  ice.  No 
functional  relationship  between  the  film  density  and  the  act¬ 
ual  surface  temperature  is  known,  since  the  IR  scanner  was  not 
calibrated.  All  that  can  be  said  is  that  the  range  of  thin 
ice  to  multi-year  ice  results  in  a  monotonic  increase  in  film 
density  from  approximately  1.0  to  2.0  density  units  for  all  five 
scans.  Since  this  range  is  typically  in  the  linear  portion  of 
the  D-log  E  or  Hur ter-Dri f f ield  curve  (for  example,  see  Thompson, 
1966,  p.  248),  it  can  be  assumed  that  the  density/ice  thick¬ 
ness  relationship  is  at  least  constant  among  the  five  scans,  and 
this  is  sufficient  to  evaluate  variations  of  skewness  of  the 
density  (or  ice  thickness)  distribution  for  these  five  samples. 

5.3  Correlation  of  Under-Ice  Roughness  with  Skewness  of 

Surface  Temperature  Distribution _ 

Figures  17-21  show  positive  images  of  the  negative 
imagery  actually  scanned  for  the  present  analysis.  Beneath 
each  image  is  a  histogram  of  the  film  density  values  derived 
from  the  scan.  Where  there  is  multi-year  ice  only  in  the 
imagery,  the  histogram  of  density  values  suggests  a  single  nor¬ 
mal  distribution.  Overlapping  normal  distributions  can  be  seen 
when  first-year  or  thin  ice  is  observed  in  the  imagery. 

For  this  work,  only  the  statistics  of  the  multi-year 
ice  distribution  are  used.  Where  there  is  an  overlapping  first 
year  distribution,  the  technique  of  subtracting  the  estimated 
first-year  ice  histogram  from  the  multi-year  histogram  is  used 
(see  for  example,  Mack,  1967,  p.  140).  From  the  histogram  of 
the  multi-year  ice,  the  skewness,  s,  or  the  amount  that  the 
'bell-shaped"  curve  leans  to  the  left  or  the  right,  was  determined 
as  follows: 

Courtesy  Research  institute,  U.S.  Army  Engineer  Topographic 
Laboratories,  Fort  Belvoir,  Virginia. 
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gram  of  film  density  values  is  displayed  below.  The  two  minor  distributions  to  the  left 
sociated  with  thin  ice  and  first-year  ice.  The  major  distribution  is  associated  with  mu 
ice  and  has  a  skewness  of  -0,830. 
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fl  +  f2  + 


where 


=  the  r  th  moment  about  the  mean 
S  =  skewness 


(  5) 


=  the  film  density  values  (the  center  of  each  class 
interval  is  chosen) 

f^  =  the  frequency  of  occurrence  of 

The  values  of  the  skewness  as  computed  from  (4)  and  (5) 
and  from  the  histograms  shown  in  Figures  17-21  are  listed  in 
Table  3.  To  be  consistent  with  the  two  examples  cited  earlier, 
we  are  defining  positive  skewness  as:  a  skewing  of  the  tempera¬ 
ture  distribution  toward  higher  temperatures,  and  negative  skew¬ 
ness  as  a  skewing  toward  lower  temperatures. 


When  the  values  of  skewness  of  the  distribution  of  ice 
surface  temperatures  are  compared  with  the  measured  under-ice 
roughness,  i.e,,  standard  deviation  of  ice  depth,  once  again  a 
clear  functional  relationship  exists,  as  shown  in  Figure  22. 

The  distribution  of  surface  temperature  measurements  of  sea  ice 
skews  toward  the  colder,  or  thicker  ice  as  deformation  increases. 
In  the  range  of  standard  deviation  of  ice  depth  covered  in  this 
example,  i.e.,  between  4.3  and  4.8  m,  the  relationship  is  linear 
and  negatixe. 


6.  Discussion  of  the  Data 


In  all  three  examples,  when  the  RMS  ice  depth  range  was  be¬ 
tween  4  and  8  m,  corresponding  to  a  standard  deviation  of  ice 
depth  ranging  between  2  and  6  m,  there  is  a  strong,  negative 
linear  correlation  between  the  skewness  of  the  temperature  dis¬ 
tributions,  whether  measured  from  satellites  or  aircraft,  and 
the  under-ice  roughness  measured  by  submarine  upward-looking 
sonar . 

Although  mathematically  a  parabola  appears  to  produce  the 
best  fit  for  all  the  data  points  used  in  the  example  presented 
in  Section  4,  the  physical  significance  of  a  parabolic  correla¬ 
tion  when  the  RMS  ice  depth  is  less  than  4  m  (3.5  m  mean  ice 
depth)  is  as  yet  unknown.  Above  4  m  RMS  ice  depth,  or  2  m 
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TABLE  3 


AREA 
( from  ] 


Comparison  of  standard  deviation  about  the  mean  ice 
depths  for  the  10  97-km  segments  shown  in  Figure  15 
with  the  skewness  of  the  film  density  distribution 
shown  in  Figures  17-21. 


Figure  15) 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


STANDARD  DEVIATION 
(ice,  in  m) 

3.995 

4.266 

4.557 

4.345 

4.783 

4.677 

4.617 

4.554 

4.057 

4.735 


SKEWNESS 
(TIR  histograms) 


-0.475 
-0.983 
-0.830 
-0.  713 
-0.580 
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SKEWNESS  OF  IR  DATA  DISTRIBUTION 


standard  deviation  (representing  most  of  the  deformed  ice  in  the 
Arctic),  the  correlation  is  simply  explained:  with  increasing 
ridging,  the  colder  the  ice  appears  relative  to  the  statistical 
mode  of  the  ice  temperature  distribution.  This  portion  of  the 
cjrrelation  curve  is  thus  consistent  with  the  other  analyses 
shown  in  Figures  8  and  22.  In  fact,  from  4  m  RMS  ice  depth  and 
deeper,  a  linear  correlation  as  was  used  in  Figures  8  and  22  can 
equally  well  fit  the  data.  It  is  highly  significant  that  a 
functional  relationship  exists  between  the  satellite  or  airborne 
TIR  and  the  under-ice  roughness  data. 

As  can  be  seen  from  Figure  2,  most  of  the  SSN  GURNARD  ice 
depth  data  falls  between  3.5  and  4.5  m  RMS.  This  is  evident  in 
Figure  14  where  most  of  the  useable  data  falls  in  the  region 
where  the  slope  of  the  correlation  curve  is  close  to  zero.  In 
this  region,  examined  with  the  other  data  excluded,  it  is  un¬ 
certain  whether  there  is  real  correlation,  i.e.,  is  the  para¬ 
bolic  function  correct,  or  is  there,  in  fact,  no  correlation. 

Even  with  the  uncertainty  of  the  physical  significance  of 
the  data  in  the  3. 5-4.0  m  RMS  ice  depth  region  it  is  clear  that 
statistically,  the  chances  are  excellent  that  there  is  a  rela¬ 
tionship  between  the  skewness  of  the  satellite  TIR  data  and  the 
RMS  under-ice  roughness  data.  Such  directional  correlations, 
consistent  with  the  actual  ice  drift,  are  highly  unlikely  to 
occur  by  chance.  In  view  of  all  of  the  above,  it  must  be  con¬ 
cluded  that  the  first  set  of  correlations  observed  with  the 
SSN  SARGO  data  taken  in  1960  and  satellite  data  taken  in  1974 
were  not  an  accident;  there  appears  to  be  a  real  correlation 
that  has  important  ramifications  for  predicting  under  ice 
roughness  for  acoustic  propagation  in  the  Arctic. 

7.  Conclusions 


Although  different  units  of  measurement,  as  dictated  by 
the  three  available  data  sets,  were  used  in  the  correlations 
to  date,  one  thing  is  clear  from  the  data  presented  here.  For 
under  ice  roughness  measured  in  terms  of  RMS  ice  depth  between 
the  ranges  of  4  m  and  8  m,  or  in  terms  of  standard  deviation 
about  the  mean  in  ranges  of  2  m  and  6  m,  there  is  a  good 
negative  linear  correlation  with  the  skewness  of  the  temperature 
distribution  of  the  corresponding  ice  surface,  where  increasing 
negative  skewness  is  defined  as  the  temperature  measurements 
skewing  toward  lower  temperatures. 

Experimental  evidence  reveals  a  definite  relationship.  To 
define  the  relationship  numerically,  the  parameters  being  cor¬ 
related  must  be  expressed  consistently.  This  is  easily  done 
with  the  under-ice  parameter  which  is  measured  in  a  straight¬ 
forward  manner,  and  can  be  expressed  as  standard  deviation  or 
RMS  ice  depth,  and  appears  conservative  over  time  (LeSchack 
and  Chang  1977). 


On  the  other  hand,  the  skewness  of  the  surface  temperature 
distribution,  which  appears  to  be  the  most  sensitive  measure¬ 
ment  of  variations  in  ice  surface  roughness  found,  has  a  numeri¬ 
cal  range  that  varies  with  the  number  of  measurements  used  in 
its  determination  and  most  probably,  the  time  of  year  the  meas¬ 
urements  were  made.  As  a  result,  it  is  likely  that  a  given 
under-ice  roughness  predictive  equation  must  be  developed  from 
T1R  data  recorded  at  essentially  the  same  time. 

It  is  concluded  that  the  under  ice  roughness  of  multi¬ 
year  ice  may  be  predicted  from  airborne  or  satellite  TIR  measure 
ments.  An  universal  predictive  equation  has  yet  to  be  developed 
however,  since  sufficient  data  are  not  presently  available. 
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APPENDIX  A 


Statistical  Analysis  of  the  VHRR  TIR  Satellite  Data 

The  satellite  TIR  data  is  broken  up  into  basic  32  x  32  km 
arrays  as  seen  in  Figure  6.  Each  32  x  32  array  is  composed  of 
1024  measurements  of  radiant  temperature  (in  the  16  x  16  km 
arrays  used  in  Section  4  of  the  report,  256  measurements  are 
contained).  Like  many  measurements,  the  temperature  measure¬ 
ments  of  multi-year  ice  appear  to  fall  into  a  Gaussian,  or 
normal  distribution.  Mathematically  we  can  express  deviations 
from  the  familiar  "bell-shaped”  curve  by  skewness  (whether  the 
bell  leans  to  the  left  or  the  right)  and  by  kurtosis,  the  peaked¬ 
ness  of  the  curve.  These  terms  are  expressed  as  follows: 


S  =  M3/(M23)l/2 
k  =  M„/M22 

and 

f i ( X i -M ) r  +  f  2 (X2-M)r  +...+  fN(XN-M)r 

M  =  - - - 

r  fi  +  fa  +...+  fN 

where 

M  =  the  r  th  moment  about  the  mean;  M  =  the  mean 
r 

S  =  skewness 
K  =  kurtosis 

X^  =  the  radiant  temperature  values 
=  the  frequency  of  occurrence  of  X^ 


N  =  the  total  number  of  values  per  matrix 
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is.  abstract  This  reportj  based  on  empirical  data,  concludes  that  a  correlation  has  been 
found  between  easily  obtainable  sea  ice  surface  temperature  and  under-ice  roughness 
data  which  are  obtainable  only  at  great  expense.  Under-ice  roughness  is  valuable  in 
evaluating  acoustic  attenuation  beneath  the  Arctic  ice  and  is  expressed  in  terms  of 
either  root-mean-square  (RMS)  ice  depth  or  standard  deviation  about  the  mean  ice  depth# 
both  of  which  are  closely  correlated.  By  showing  a  functional  relationship  between  the 
skewness  of  the  surface  temperature  distribution  as  derived  from  NOAA  VHRR  Satellite 
thermal  infrared  data  and  the  under-ice  roughness,  the  way  appears  clear  to  make  a 
chart  of  Arctic  under-ice  roughness  for  Arctic  acoustic  programs  and  for  nuclear  sub- 
ir  diners . 

. ^  Specifically,  a  comparison  of  the  skewness  of  the  temperature  distribution  de¬ 
rived  from  NOAA  VHRR  thermal  infrared  satellite  data,  of  the  Beaufort  Sea  taken  in  April 
1974  was  made  with  the  RMS  under-ice  roughness  derived  from  under-ice  profile  data  re¬ 
corded  by  the  SSN  SARGO  in  February  1960  taken  at  essentially  the  same  locations.  Five 
data  sets  were  chosen  ranging  from  the  highly  deformed  and  hence  thicker  ice  area  near 
the  Canadian  Archipelago  to  thinner  ice  off  the  coast  of  Alaska  near  Prudhoe  Bay.  A 
very  high  negative  correlation  was  observed. 

In  a  second  study,  under  ice  data  recorded  in  April  1976  by  the  SSN  GURNARD 
was  correlated  with  the  skewness  of  temperature  distributions  derived  from  NOAA  VHRR  IR 

data  recorded  in  March  1976  over  nominally  the  same  area  of  the  Beaufort  Sea. _ The  SSN 

GURNARD  cruise  track  of  some  780  N.  miles  was  divided  into  63,  10  N.  mile  (16  km)  sec-^ 
tions.  The  NOAA  VHRR  IR  data  were  similarly  grouped  into  10  N.  mile  squares  that 
nominally  corresponded  to  the  positions  from  where  RMS  under  ice  data  were  gathered. 
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ABSTRACT  continued 

Twenty  useable  skewness  values  were  obtained  with  the  RMS  under  ice  data.  Because  of 
the  relatively  high  resolution,  i.e. ,  10  N.  miles,  some  data  set  misregistration  was 
anticipated,  but  by  moving  the  satellite  data  matrix  about  10  NM  to  the  west,  north-  J 
west  or  southwest,  a  good  registration  was  found  and  a  good  statistical  correlation 
was  obtained.  This  slight  matrix  shifting  appears  significant  because  it  is  con¬ 
sistent  in  both  direction  and  magnitude  with  the  actual  measured  ice  drift  between 
March  and  April  1976.  Also,  shifting  in  any  of  the  other  three  quadrants  shows  no 
correlation . 

x  A  third  study  was  then  conducted  that  shows,  perhaps  more  graphically  than 
the  others,  the  correlation  of  under- ice  data  recorded  by  the  British  nuclear  sub¬ 
marine  HMS  SOVEREIGN  between  18-21  October  1976  with  airborne  IR  data  recorded  dur-  ( 
ing  the  same  period  over  the  submarine  track  by  a  Canadian  Forces  Argus  aircraft. 

The  airborne  IR  data  are  recorded  in  the  form  of  imagery  on  a  scale  large  enough'  to 
see  the  individual  features  of  the  ice  surface.  A  scanning  microdensitometer  gen¬ 
erated  histograms  of  film  density  values  that  are  proportional  to  ice  surface  temp¬ 
eratures.  The  skewness  of  the  histograms  were  correlated  with  the  standard  devia¬ 
tions  about  the  mean  ice  depth  for  the  associated  submarine  track. 

— ‘  » 

In  all  three  examples,  when  the  RMS  ice  depth  range  was  between  4  and  8m,  ( 

corresponding  to  a  standard  deviation  of  ice  depth  ranging  between  2  and  6  m,  there 
is  a  strong,  negative  linear  correlation  between  the  skewness  of  the  temperature  dis-  <• 
tributions,  whether  measured  from  satellites  or  aircraft,  and  the  under-ice  rough-  (| 
ness  measured  by  submarine  upward- look  ing  sonar.  jj 
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